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Abstract. We give a constructive and exhaustive definition of Kochen-Specker (KS) 
vectors in a Hilbert Space of any dimension as well as of all the remaining vectors 
of the space. KS vectors are elements of any set of orthonormal states, i.e., vectors 
in ri-dim Hilbert space, Ti. n , n > 3 to which it is impossible to assign ls and 0s in 
such a way that no two mutually orthogonal vectors from the set are both assigned 
1 and that not all mutually orthogonal vectors are assigned 0. Our constructive 
definition of such KS vectors is based on algorithms that generate MMP diagrams 
corresponding to blocks of orthogonal vectors in K n , on algorithms that Single out 
those diagrams on which algebraic 0-1 states cannot bc defined, and on algorithms 
that solve nonlinear equations describing the orthogonalities of the vectors by means 
of statistically polynomially complex interval analysis and self-teaching programs. The 
algorithms are limited neither by the number of dimensions nor by the number of 
vectors. To demonstrate the power of the algorithms, all 4-dim KS vector Systems 
containing up to 24 vectors were generated and described, all 3-dim vector Systems 
containing up to 30 vectors were scanned, and several general properties of KS vectors 
were found. 
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1. Introduction 

Recently proposed experimental tests of the Kochen-Specker (KS) theorem [H Ej, 
skepticism on the feasibility of such experiments [3-7], positive experiments recently 
carried out [Hj, and recent theoretical elaborations on the theorem [9-22] prompted a 
renewed interest in the KS theorem. 

The KS theorem proves that there is a set of measurements that can be carried 
out on a finite dimensional quantum system in such a way that if one assumed that 
the values of measured observables are completely independent of all other observables 
that can be measured on the same System, then one would run into a contradiction. 
Hence, a quantum system cannot posses a defmite value of a measurable property prior 
to measurement, and quantum measurements (essentially detector clicks) carried out on 
quantum Systems cannot be ascribed predetermined values (say and 1). To arrive at 
the claim, one considers an orthonormal set of states {ipi, . . . , ip n }, i.e., vectors in n-dim 
Hilbert space, TC n , n > 3. Projectors onto these states satisfy: Y17=i^ = ^ where 
Pi = ipiipj. Now, Kochen and Specker proved [2H| that there is no function / : TC — > R 
satisfying the Sum Rule Yli=i fi-^d = f(Yli=i Pi) = /CO f° r a ^ se ts °f projectors Pj. 
Hence, there is at least one set of projectors {Pj, P/, . . .} and the corresponding set of 
vectors {tßi,^, ■ ■ .} for which the Sum Rule is not satisfied. Choosing /(Pj) G {0, 1} 
(/(0 = 1); the theorem amounts to the following claim: In 7i n , n > 3, it is impossible 
to assign ls and 0s to all vectors from such a set — which we call a KS set — in such a 
way that [24J: 

(i) No two orthogonal vectors are both assigned the value 1; 

(ii) In any subset of n mutually orthogonal vectors, not all of the vectors are assigned 
the value 0. 

All the vectors from a KS set, as defined above, we call KS vectors. KS vectors in 
each KS set form subsets of n mutually orthogonal vectors. We arrive at one subset from 
another by a series of rotation in 2-dim planes around (n— 2)-dim subspaces as explained 
in See. Thus, any two subsets share at least one vector which is orthogonal to all other 
vectors in both subsets and in an n-dim space, two subsets can share up to n— 2 vectors. 
The KS vectors correspond to the directions of the quantisation axes of the measured 
eigenstates within experiments which have no classical counterparts, and when we speak 
of finding KS vectors we mean finding these directions. We stress here that it is not our 
aim to give yet another proof of the KS theorem but to determine the class of all KS 
vectors from an arbitrary H, n as well as the class of all non-KS vectors, i.e., vectors from 
the remaining sets of vectors from HJ 1 . By the class of non-KS vectors we mean vectors 
that allow 0-1 states and that correspond to the directions of the quantisation axes of 
the measured eigenstates within experiments which do have classical counterparts and 
when we speak of finding non-KS vectors we mean finding the latter directions. 

The original KS theorem [23J made use of 192 (claimed 117) 3-dim vectors. 
Subsequent attempts to reduce the number of vectors gave the following minimal results 
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(usually called records): Bub's System contains 49 vectors (claimed 33) [23] > Conway- 
Kochen's has 51 (claimed 31) [2H1 P- 114], and Peres' System has 57 (claimed 33) [27] 
3-dim vectors f ; Kernaghan's System contains 20 4-dim vectors with the smallest loops 
(see the defmition below) of size two [2H]; Cabello's System has 18 4-dim vectors with 
the smallest loops of size three (301 , etc. Reducing the number of vectors is important 
for devising experimental Setups [Tüj . especially so as recently a Single qubit KS scheme 
was formulated |9, by means of auxiliary quantum Systems (ancillas) of the measuring 
apparatus and subsequently connected with the original KS formulation [TQj . On the 
other hand, knowing the class of all KS vectors is important for better theoretical insight 
into quantum theory and possibly designing quantum Computers. However, no general 
method for constructing sets of KS vectors has been proposed so far and the aim of this 
paper is to give one. In doing so we will follow the ideas put forward in [3T | 132 ] 133]. 

So far, KS vectors have been constructed either by means of partial Boolean algebras 
and orthomodular lattices [23 EU EHl EU] , by direct experimental proposals [U 121 HE], or 
by combining rays in R n |23 | l2T | I2TJ]. These approaches have two disadvantages: first, 
they depend on human ingenuity to find ever new examples and "records," and second, 
their complexity grows exponentially with increasing numbers of dimensions and vectors. 
For example, lattices of orthogonal n-tuples have 2 n elements (Hasse diagrams) jlü] and, 
on the other hand, the complexity of nonlinear equations describing combinations of 
orthogonalities also grows exponentially. 

As opposed to this, we are able to give algorithms for generation of all the 
equations that have KS vectors as their Solutions and to effectively solve them (up 
to a reasonably chosen number of vectors and dimensions — limited only by the speed of 
today's Computers) in a way that is essentially of a statistically polynomial complexity. 
We first recognise that a description of a discrete observable measurement (e.g., spin) 
in Ti. n can be rendered as a 0-1 measurement of the corresponding projector along the 
vector in ¥L n onto which the projector projects. Hence, we deal with orthogonal triples 
in M. 3 , quadruples in M. 4 , etc., which correspond to possible experimental designs, and 
to find KS vectors means finding such n-tuples in M. n . 

The orthogonalities of vectors within these n-tuples can be described by nonlinear 
equations of type given in Eq. Q that have Solutions. There are however billions of 
such nonlinear equations that have no Solutions even for the smallest KS sets. And 
their number grows exponentially with the increase of both the number of KS vectors 
and the dimension of their space. So, we established a one-to-one correspondence 
between nonlinear equations and graphs (MMP diagrams). We can handle graphs 
exponentially faster than nonlinear equations but there are nevertheless billions of them. 
Therefore, we designed a self-teaching generation algorithm for MMP diagrams: graphs 
containing subgraphs that correspond to equations that cannot have a Solution are not 
generated. This reduces the generation complexity to a statistically polynomial one 

f The reasons why Kochen-Specker's, Bub's, Conway-Kochen's, and Peres' Systems should be 
considered as 192, 49, 51, and 57 and not as 117, 33, 31, and 33 vector Systems, respectively, are 
given in See. ß3-(xi) in aecordance with the results independently obtained by J.-Ä. Larsson. [2H| 
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and the time required for obtaining the MMP diagrams corresponding to Systems of 
nonlinear equations with Solutions from billions of years to hours and days. 

To switch back from the MMP diagrams to nonlinear equations to solve them at 
this stage would again take quite some time. Therefore we defined the notion of an 
algebraic dispersion-free state (0-1 state) on MMP diagrams. It turns out that only 
a small percentage of the obtained MMP diagrams cannot have 0-1 states. Their 
direct verification is again of exponential complexity. So, we developed algorithms with 
backtracking that discard MMP diagrams with 0-1 states and whose complexity turns 
out to be statistically polynomial. 

The diagrams finally obtained correspond to candidate sets of nonlinear equations 
that contain KS sets provided the equations have real Solutions. Algorithms for solving 
nonlinear equations, such as Gröbner basis and homotopy, are also mostly of at least 
exponential complexity and have been tested without success on representative Systems. 
Therefore we designed new ones based on interval analysis and Ritt's characteristic set 
calculations and were able to reduce their complexity to a statistically polynomial one. 
This rounds up the constructive and exhaustive definition of KS sets and vectors and 
makes their generation feasible for reasonably chosen numbers of vectors and dimensions. 

The paper is organised as follows. In See. El MMP diagrams are defined and the 
algorithms as well as the programs for their generation are presented. In See. El we give 
the algorithm and program for finding whether MMP diagrams can be assigned a set of 
dispersion-free 0-1 states and determining the latter sets when there is at least one set 
of 0-1 states together with the smallest MMP diagrams that do not allow 0-1 states. 
In See. we establish a link between MMP diagrams that do not allow 0-1 states and 
Systems of nonlinear equations whose Solutions are the KS vectors. We then give the 
algorithms and methods for solving the equations in a statistically polynomial time. In 
See. we present the new results we obtained. 

2. MMP diagrams 

We start by describing vectors as vertices (points) and orthogonalities between them as 
edges (lines connecting vertices), thus obtaining MMP diagrams jSH ESI ES] which are 
defined as follows: 

1. Every vertex belongs to at least one edge; 

2. Every edge contains at least 3 vertices; 

3. Edges that intersect each other in n — 2 vertices contain at least n vertices; 

Isomorphism-free generation of MMP diagrams follows the general principles 
established by [36J, which we now recount briefly. 

Deleting an edge from an MMP diagram, together with any vertices that lie only 
on that edge, yields another MPP diagram (perhaps the vaeuous one with no vertices). 
Consequently, every MMP diagram can be construeted by starting with the vaeuous 
diagram and adding one edge at a time, at each stage having an MMP diagram. 
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We can represent this process as a rooted tree whose vertices correspond to MMP 
diagrams whose vertices and edges have unique labels. The vacuous diagram is at the 
root of the tree, and for any other diagram its parent node is the diagram formed by 
deleting the edge with the highest label. The isomorph rejection problem is to prune 
this tree until it contains just one representative of each isomorphism class of diagram. 
This can be achieved by the application of two mies. 

Given a diagram D, we can identify the valid positions to add a new edge such that 
Conditions 3-4 are enforced. According to the symmetries of D, some of these positions 
are equivalent. The first rule is that exactly one position in each equivalence class of 
positions is used; a node in the tree formed by adding an edge in any other position is 
deleted together with all its descendants. 

To understand the second rule, consider a diagram D' with at least one edge. We 
label the edges of D' in a canonical order, which is an order independent of any previous 
labelling. Then we dehne the major class of edges as those that are equivalent under 
the symmetries of D' to the edge that is last in canonical order. The second rule is: 
when D' is constructed by adding an edge e to a smaller diagram, delete D' (and all its 
descendants) unless e is in the major class of edges of D' . 

According to the theory in [3ü]> application of both rules together is sufficient: 
exactly one diagram from each isomorphism class remains in the tree. Our implementa- 
tion used nauty |3Zj for Computing symmetries and canonical orderings. The method 
allows for very efficient parallelisation of the computation. A generation tree for MMP 
diagrams with 9 vertices and the smallest loop of size 5 is shown in the Fig. ^ 

1 2 3 




Figure 1. An example of a generation tree for connected MMP diagrams: 9 vertices 
and the smallest loop of size 5 (for 9 vertices a loop cannot be formed; the first loop 
appearswith 10 vertices: 123, 345, 567, 789, 9A1). Cf. [331 EU 

MMP diagrams with three vertices per edge and with smallest loops (edge polygons) 
of size five graphically resemble Greechie diagrams [ÜB] • Greechie diagrams are a handy 
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way to draw Hasse diagrams that represent orthomodular lattices. The complexity of 
Hasse diagrams grows exponentially with increasing dimensions and the smallest loops 
of the corresponding are of size 5, while MMP diagrams allow loops of size 4, 3, and 2 
(2 edges share at least 2 vertices). Besides, it would be quite a challenge to find a direct 
lattice representation of KS vectors. 

We denote vertices of MMP diagrams by l,2,..,A,B,..a,b,.. By the above 
algorithm we generate MMP diagrams with chosen numbers of vertices and edges and 
a chosen minimal loop size. E.g., in the examples given in Fig. Elwe generate diagrams 
with 4 vertices within an edge and minimal loops of size 2 and 3. Our programs handle 
diagrams with up to 90 vertices, but this limit could easily be extended. 

3. Algebraic states on MMP diagrams 

To find diagrams that cannot be ascribed 0-1 values we apply an algorithm which we call 
statesOl. The algorithm is an exhaustive search of MMP diagrams with backtracking. 
The criterion for assigning 0-1 (dispersion-free) states is that each edge must contain 
exactly one Vertex assigned to 1, with the others assigned to 0. As soon as a vertex 
on an edge is assigned a 1, all other vertices on that edge become constrained to 0, 
and so on. The algorithm scans the vertices in some order, trying then 1, skipping 
vertices constrained by an earlier assignment. When no assignment becomes possible, 
the algorithm backtracks until all possible assignments are exhausted (no Solution) or 
a valid assignment is found. In principle the algorithm is exponential, but because the 
diagrams of interest are tightly coupled, constraints build up quickly. For the ränge of 
diagram size in our study, we found that the average time per diagram appeared to grow 
polynomially with the diagram size. 

To implement the algorithm we wrote a program that selects MMP diagrams with 
3 and 4 vertices per edge on which 0-1 states cannot be defined. The smallest such 
diagrams are given in Fig. El 




Figure 2. Smallest MMP diagrams without 0-1 states: (1) 4 vertices per edge: (a) 
loops of size 2: 6 vertices — 3 edges; (b) loops of size 3: 10-5; (c) loops of size 4: 22-11; 
(2) 3 vertices per edge: (d) loops of size 5: one of two 19-13; the other is shown in 
Fig. 2(b) of [33!- 
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• 3 vertices per edge 

7 vertices — 5 edges (smallest loops of size 3): 123,345,561,275,476 (triangle); 
15-11 (4): 123,345,567,789,9AB,BC1,CD6,2DA,2E8,4FA,CEF (hexagon), 
123 , 345 , 567 , 789 , 9 AB , BCD , DE1 , 4AE , 28C , 2FA , 6FD (heptagon) ; 
19-13 (5): 123,345,567,789,9AB,BCD,DEF,FG1,2IA,6IE,4HC,8JG,HIJ Fig.0(d), 
123 , 345 , 567 , 789 , 9 AB , BCD , DE1 , EI7 , 2F9 , 4GB , I JG , F JH , CH6 (heptagon) ; 

• 4 vertices per edge 

6-3 (smallest loops of size 2): 1234,2356,1456 Fig.0(a); 

10-5 (smallest loops of size 3): 1234, 4567, 7891, 35A8,29A6 Fig.H(b); 

22-11 (smallest loops of size 4): 

1234 , 4567 , 789A , ABCD , DEFG , GHI 1 , FJK5 , H JMC , 3KL8 , IBL6 , 29ME . Fig. H(c) , 
1234 , 4567 , 789A , ABCD , DEF1 , FGH5 , EMJ6 , 2GLC , 3IJ8 , HIKB , MLK9 (pentagon) ; 

38-19 (5) : 1234 , 1567 , 289A , 5BCD , 8BEF , 3GHI , 6JKL, GJMN , CHOP , EMQR , OQST , RUVW , 



Further details of the algorithm and the program statesOl will be given elsewhere. jl2] 
4. Kochen-Specker vectors 

To find KS vectors we follow the idea put forward in [SU ESI an d proceed so as to 
require that their number, i.e. the number of vertices within edges, corresponds to the 
dimension of ~R n and that edges correspond to n(n — l)/2 equations resulting from inner 
products of vectors being equal to zero which means orthogonality. So, e.g., an edge of 
length 4, BCDE, represents the following 6 equations: 



2L B 


• ac 


= (lBlO>Cl + Oß20C2 + ClB3 a C3 + ÖB4ÖC4 


= 0, 


& B 


■ &d 


= ÜBlClDl + OB2QD2 + Q_B3ÖD3 + Oß4ÖD4 


= o, 


B. B 


• a^ 


= dBX&EX + aß2 a E2 + Oß3 a £;3 + ÖB4Ö.E4 


= 0, 


a c 


• a^ 


= üciCLDl + CtC2 a D2 + a C3 a D3 + &C± a DA 


= 0, 


a c 


• a^ 


— O-ClO-El + «C2Ö_E2 + OC3 a £3 + 0>Ci a EA : 


= 0, 


an 


■ a E 


= aoid-El + 0>D20'E2 + a D3 a B3 + Q-DA^EA 


= 0. 





(1) 



Each possible combination of edges for a chosen number of vertices corresponds to a 
System of such nonlinear equations. A Solution to Systems which correspond to MMP 
diagrams without 0-1 states is a set of components of KS vectors we want to find. Thus 
the main clue to finding all KS vectors is the exhaustive generation of all MMP diagrams 
as given in See. then picking out all those diagrams that cannot have 0-1 states as 
presented in See. El establishing the correspondence between the latter diagrams and 
the equations for the vectors as shown in Eq. and finally solving the Systems of the 
so obtained equations. 



4UXY , 9SZa , FIbc , KTXb , 7VZc , ALPW , DNYa (dodecagon) . 




Kochen- Specker Vectors 



8 



In practice, we actually merge these four stages so as to avoid generating those 
diagrams that cannot have a Solution, (j For Systems of equations of type given by 
Eq. (JH) that do have Solutions that do not allow 0-1 states, such Solutions are KS 
vectors that correspond to vertices of MMP diagrams. Mutually orthogonal vectors 
correspond to edges, and connected edges, i.e., MMP diagrams themselves correspond 
to the Systems of equations. For instance, in the connected edges 1234,4567 vectors 
1,2,3,4 and 4,5,6,7, are mutually orthogonal and 4567 is obtained from 1234 by 
4-dim rotations (1234 and 4567 are connected by 4). A general 2-dim rotation is a 
rotation by an angle around a fixed point in the same plane. A general 3-dim rotation is 
a rotation in a 2-dim plane by an angle around a fixed axis perpendicular to this plane. 
So, we define a general 4-dim rotation as a rotation in a 2-dim plane by an angle around 
a fixed 2-dim plane. What is common to all these rotations is that they always take 
place in a 2-dim plane. Hence, we define an n-dim rotation as a rotation in a 2-dim 
plane by an angle around a fixed (n — 2)-dim subspace [H]. This also explains the case 
of a smallest loop of size 2 in the 4-dim case. E.g., we arrive at 4561 from 1234 by a 
rotation in the 2-dim plane determined by the vectors 2 , 3 (and also by the vectors 5 , 6) 
around the plane determined by the vectors 1 , 4. 

Finding KS vectors is not a well-posed problem in terms of solving, though. Indeed 
if V is a KS vector then AV is also a KS vector for any non-zero scalar A. Furthermore, 
if S is a set of KS vectors, then TZS is also such a set for any arbitrary rotation matrix 
TZ. We may simplify the problem by considering only unit vectors (i.e. vectors whose 
Euclidean norm is 1 and hence vectors whose components have a value in the ränge [- 
1,1]). To avoid the rotation problem, we may assume that one n-tuple is the orthonormal 
basis of R n . Under these assumptions, some of the orthogonality equations simplify. 
E.g., if 1234 is the basis of M 4 with 1=[0, 0,0,1] then 1567 indicates that the fourth 
components of 5,6,7 are 0. The non-collinearity constraints also plays an important 
role. E.g., 1235 is not a possible n-tuple as three components of 5 would be and hence 
5 would be collinear with 4. 

This has prompted us to develop a preliminary pass, which allows elimination of 
n-tuples that cannot lead to a Solution. Consider a System of m 4-dim vectors. The 
preliminary pass makes use of an m x 4 table T, called the 0-table, with an entry set 
to 1 when a vector component cannot be 0. For example, if vector j has components 
[0^1,0,0,0^4], then neither aji nor 0,^4 can be (otherwise the vector will be collinear 
with one of the vectors of the basis) and T[j, 1] = T[j, 4] = 1. The preliminary pass 
selects a set of four 4-dim vectors as the basis of 1R 4 . It then applies a set of simplification 
rules on the the orthogonality equations. For example, if the equation is ajk Ojfc = and 

k] — 1 then a^, is set to 0. Each time a vector component value is determined the 
0-table is updated and the preliminary pass is restarted. The process will stop when no 
further simplification may be performed or when a constraint violation occurs (e.g., one 
equation implies that should be while the 0-table indicates that this component 

jj This merging is crucial. Without it we would not be able to reduce the exponential complexity of 
the problem to the statistically polynomial one. 
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cannot be 0) in which case the System cannot have a Solution. The simplification rules 
used in the preliminary pass depend on the space dimension. 

The preliminary pass has been implemented as a C program that has been added 
as a filter in the generation program. For avoiding the exponential growth of the 
number of generated MMP diagrams it is essential that the candidate KS-sets should be 
generated incrementally i.e. that the program generates sequentially all Systems starting 
with a given m n-tuples before modifying the mth n-tuple. By using this incremental 
generation during the preliminary pass determines that an initial set of m n-tuples has 
no Solution and that no further Systems starting with this set will be generated. E.g., 
for 18 vectors and 12 quadruples, without such a filter we would generate > 2.9 ■ 10 16 
Systems — what would require more than 30 million years on a 2 GHz CPU — while the 
filter reduces the generation to 100220 Systems (obtainable within < 30 mins on a 2 
GHz CPU). Thereafter statesOl gives us 26800 Systems without 0-1 states in < 5 secs. 

For the remaining Systems, two solvers have been developed. One is based on 
a specific implementation of Ritt characteristic set calculation [13]. Assume that a 
vector V=[ay 1; 0, 0, ay 4 ] (this implies that a vi ,av4: cannot both be 0) is orthogonal 
to W=[awi, 0, aw2, a>w4\- From the orthogonality condition awi a vi + o-Wi oy4 = 0, 
we deduce that awi = —awi avi/ a vi (as a V i cannot be 0) and that a WA cannot be 
0. This Information is propagated to the other equations and, as for the preliminary 
pass, simplification rules are applied to the equations, allowing us to determine 
further unknowns and to update the 0-table. The process is repeated until no 
further unknowns can be determined or until a constraint violation occurs. If no 
violation occurs, we will usually get a set of remaining equations that is quite simple 
and that allows us to determine all Solutions. This solver has been implemented 
using the symbolic computation Software Maple. E.g., for system (a) in Fig. El 
1234 , 4567 , 789A , ABCD , DEFG , GHI 1 , 35CE , 29BI , 68FH, we get (in < 10 secs on a 2 GHz 
CPU) the remaining set of 10 equations: 2üg 2 = 2a 2 ci = 2üq 3 = 2a\± = ^a 2 EA = 1, 
2o| 3 = 2a 2 I2 = 2ag 4 = 4a^ 2 = 2aj 1 = 1 from the roots of which we may deduce the 
other vector components (e.g., we get 6 as a 62 [0, 1, —2054094053, ag 4 ]). The drawback 
of this approach is that it does not always allow us to completely solve the equational 
Systems: we may end up with a system with fewer equations for which no further 
constraints can be propagated. 

Our second solver is based on interval analysis. An interval evaluation of an 
equation f(xx,..., x m ) = is a ränge F = [a, b] such that if all the unknowns Xi, . . . ,x m 
are restricted to lie within given ranges, then whatever is the values of the unknowns 
in their ränge we have a < f(x±, . . . ,x n ) < b. A simple way to calculate an interval 
evaluation is to use interval arithmetic that simply replaces all mathematical Operators 
by an interval equivalent. E.g., the interval evaluation of the orthogonality condition 
XiVi + with Xi,x 2 G [0.5,1] and yi G [0.1,0.2], y 2 G [0.2,1] is calculated as 
[0.5,1] [0.1, 0.2] + [0.5,1] [0.2,1] = [0.05,0.2] + [0.01,1] = [0.06,1.2]. Note that if the 
interval evaluation of an equation does not include then there is no value of the 
unknowns in their ränge that can cancel the equation. A box will be a set of ranges, one 
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Figure 3. Smallest 4-dim KS Systems with: (1) loops of size 3: (a) 18-9 (isomorphic to 
Cabello et al. (SOj); (b) 24(22)-13 not containing System (a), with values ^ { — 1,0,1}; 
(2) loops of size 2: (c) 19(18)-10. 

for each unknowns. 

Solving a KS-system is an appropriate problem for interval analysis, since all the 
unknowns are in the ränge [-1,1]. The set of these unknowns is the box B . The System of 
equations to be solved consists of the equations derived from the orthogonality conditions 
between the vectors and of the unitary equations that describe that each vector is a unit 
vector. 

A basic solver uses a list of boxes £ that initially has element £>o- At step 
i, the algorithm processes box Bi of C and calculates the interval evaluation of the 
orthogonality and unitary equations: if the interval evaluation of one of these equations 
does not include 0, then the algorithm will process the next box in the list. Otherwise 
two new boxes will be generated from Bi by bisecting the ränge of the box. These 
boxes will be added to the list, and the next box in the list will be processed. The 
algorithm will stop either when all the boxes have been processed (meaning the System 
has no Solution) or when the width of all the ranges in a box is less than a small value 
while the interval evaluations of all the equations still include (meaning a Solution is 
obtained). Note that the method is mostly sensitive to the number of unknowns (which 
explains why the vectors that appear only once in the KS-system should be eliminated) 
and not so much on the number of equations. On the contrary, additional equations may 
even reduce the computation time. For example, consider a triplet XiXjXk in 3D: using 
the orthogonality condition, Xk is obtained as Xk = ±X{ x Xj. We get therefore two 
possible Solutions for Xk, and additional equations will be obtained by writing that the 
Square of each component of Xk should be equal to the Square of the same component 

of X; X Xj, 

Numerous methods may be used to improve the efficiency of the basic solver 
(especially to prove that indeed a System has a Solution). We use the interval analysis 
library ALIAS f to deal emciently with the KS Systems. 

f www.inria.fr / coprin/logiciels / ALIAS / ALIAS .html 
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Interval analysis has in principle an exponential complexity, due to the bisection 
process. But it has been experimentally shown that in some cases, the practical 
complexity is only polynomial. According to our tests (over 400 billion Systems have 
been checked) it appears that the solving of the KS Systems has indeed only a statistically 
polynomial complexity. It must also be noted that the solver may be used during the 
generation of the MMP diagrams as a complement to the preliminary pass to avoid the 
exponential growth of the number of generated MMP diagrams. Indeed, for a diagram 
that has not been rejected by the preliminary pass, we may run the solver to further 
check if the diagram has a Solution. But since the solver may be relatively Computer 
intensive, we have to use an adaptive version in which the number of allowed bisections 
is limited. For example this number may be large for relatively small sub-graphs because 
determining that they don't have a Solution allows us to avoid the generation of a large 
number of diagrams. On the other hand, the number of allowed bisections will be small 
for sub-graphs whose size is close to the maximum (and consequently from which few 
diagrams will be deduced), thus avoiding increased generation time. 

We also developed a checking program that finds Solutions from assumed sets, say 
{ — 1,0, 1}, even faster (< 1 sec on a 2 GHz CPU) by precomputing all possible scalar 
products. The main algorithm scans the vertices and tries to assign unique vectors to 
them so that all vectors assigned to a given edge are orthogonal. In case of a conflict 
the algorithm backtracks, until either all possible assignments have been exhausted or a 
Solution is found. We match its exponential behaviour by scanning next those vertices 
most tightly coupled to those already scanned, helping to force conflicts to show up 
early on so that backtracking can take care of them more quickly. 

Further details of the algorithms and programs presented in this section will be 
given elsewhere. 

5. New results and conclusions 

In this paper we presented algorithms that generate and those that solve sets of 
arbitrary many Kochen-Specker (KS) vectors that are of polynomial complexity or at 
least of statistically polynomial complexity. The algorithms merge generation of MMP 
diagrams corresponding to blocks of orthogonal vectors in R n (Sec. |2J), singling out 
MMP diagrams on which 0-1 states cannot be defined (Sec. EJ), and solving nonlinear 
equations describing the orthogonalities of the vectors by means of interval analysis 
(Sec. EJ), so as to eventually generate KS vectors in a statistically polynomially complex 
way. Using the algorithms we obtained the following results: 

(i) A general feature we found to hold for all MMP diagrams without 0-1 states 
we tested is that the number of edges, b and the number of vertices that share more 
than one edge, a* satisfy the following inequality: nb > 2a*, where n is the number of 
vertices per edge. Hence, there are no KS vectors that share at least 2 of & n-tuples in 
their KS set whose number a* > In R n this means that we cannot arrive at Systems 
with more unknowns than equations when we disregard the unknowns that appear in 
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only two equations. To prove the feature for an arbitrary n remains an open problem. 

(ii) For MMP diagrams without 0-1 states with 3 vertices per edge and a < 30 as 
well as with 4 vertices per edge and a < 23 the stronger inequality holds: nb > 2a. The 
only exception to this rule we have found is the original Kochen-Specker system with 
192 vertices [see (xi) and Fig. Ej. At what a for a chosen n this inequality ceases to hold 
is an open problem. 

(iii) None of the Systems corresponding to the smallest diagrams without 0-1 states 
given in See. El and Fig. EJhas a Solution. || The smallest KS vectors that we found to 
have real Solutions are presented in Figs. El and |U 

(iv) Between the 4-dim system shown in Fig. El(a) and the one shown in 
Fig. E](a) (both with smallest loops of size 3) there are 62 Systems with loops 
of size 3, all containing the system (a), 37 of which do not have Solutions from 
{ — 1,0,1}. System (b) is the first system with loops of size 3 not containing (a): 
1234 , 4567, 789A , ABCD , DEFG , GHI 1 , FNM8 , G0L7 , HJK6 , DNK4 , AMJ1 , 35CE , B29 J. One of its 
Solutions is: 12. . . NO = {1,0,1, 1}{1,0,-2,1}{1,0,0,-1}{0,1,0,0}{0,0,1,0}{0,0,0,1}{1,0,0,0} 
{0,2,2,1} {0,2,-1,-2} {0,1,-2,2} {3,2,2,1} {1,-2,0,1} {-1,0,1,1} {1,1,0,1} {1,-1,1,0} {0,1,1,-1} 
{l,l,-l,0}{l,-l,0,-l}{l,-2,-l,0}{l,0,l,0}{0,0,l,l}{3,2,-l,-2}{l,0,-l,2}{0,2,-l,l} (which can, 
of course, easily be normalised. The system does not have a Solution from {—1, 0, 1}. 

(v) The smallest 4-dim system with the smallest loop of size 2 is the following 19-10 
one: 1234, 4567, 789A, ABCD, DEFG, GHI 1,35CE,29BI,68FH, 6781 shown in Fig. E](c). It 
contains system (a) of Fig. El and it is the only MMP system with 19 vertices which has 
a Solution from {—1, 0, 1} for the corresponding vectors. 

(vi) The two smallest 4-dim Systems with the smallest loops of size 2 
that do not contain system (a) of Fig. El are the following 20-10 ones: 
1234 , 4567 , 789A , ABCD , DEFG , GHI 1 , 68FH , 12 JI , 1J9B , 345K , 4KEC and 1234 , 4567 , 789A , 
ABCD,DEFG,GHI1,68FH,2IAK,345J,4JEC,9ABK shown in Figs. ü(a) and (b). The lat- 
ter system is isomorphic to Kernaghan's system (2H1- A Solution to the former one is: 
12. . . JK = {0,0,0,1}{1,0,0,0}{0,1,1,0}{0,1,-1,0}{1,0,0,-1}{1,1,1,1}{1,-1,-1,1}{1,1,-1,-1} 
{1,0,1,0}{0,1,0,1}{1,0,-1,0}{1,1,1,-1}{1,-1,1,1}{1,-1,-1,-1}{0,0,1,-1}{1,1,0,0}{1,-1,0,0} 
{0,0,1,0}{0,1,0,0}{1,0,0,1}. 

(vii) All 4-dim Systems with up to 22 vectors and 12 edges with the small- 
est loops of size 2 which do have Solutions from { — 1,0,1} contain at least one 
of the Systems (a) and (b) of Fig. 0] and in many cases also (a) of Fig. El 
The two smallest 4-dim Systems with the smallest loops of size 2 that contain 
neither of the latter three Systems are 22-13 Systems (c) and (d) of Fig. 0J 
1234 , 4567 , 789A , ABCD , DEFG , GHI 1 , 2ILA , 345 J , 4 JEC , 678K , 7KMG , 9 ABL , FGHM and 
1234 , 4567 , 789A , ABCD , DEFG , GHI 1 , 121 J , 345K , 678L , GML7 , 1J9B , 4KEC , FGHM. Their 

|| Still, they might be significant for other fields. E.g., the two diagrams 19-13(5) given in See. Ol [one of 
them is shown in Fig. 2 (b) of j^Hl and the other in Fig. El(d)] are equivalent to the Greechie diagrams 
with 19 atoms and 13 blocks and to our knowledge, the smallest Greechie diagram with 3 atoms per 
edge without 0-1 states known so far was the one given by Greechie [441 I45| . with 27 atoms and 18 
blocks. The System 38-19(5) from See. ßlyields the smallest Greechie diagram with 4 atoms per block. 
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Figure 4. Smallest 4-dim KS Systems with loops of size 2: (1) — not containing System 
(a) of Fig. |3| (a) 20-11; (b) 20-11 isomorphic to Kernaghan [221; (2) — containing 
ncither System (a) of Fig. |3 nor Systems (a) and (b) off this figure: (c) 22-13; (d) 
22-13. 



Solutions are: 12. . .M={1,1,0,0}{1,-1,0,0}{0,0,1,0}{0,0,0,1}{1,0,0,0}{0,1,1,0}{0, 1,-1,0} 
{1,0,0,1}{1,-1,-1,-1}{1,1,1,-1}{1,-1,1,1}{1,0,-1,0}{0,1,0,1}{1,0,1,0}{1,1,-1,-1}{1,-1,-1,1} 
{l,-l,l,-l}{0,0,l,l}{0,l,0,0}{l,0,0,-l}{l,l,-l,l}{l,l,l,l}and 12. . .M={0,0,0,1}{1,0,0,0} 
{0,1,1,0} {0,1,-1,0} {1,0,0,-1} {1,1,1,1} {1,-1,-1,1} {1,-1,1,-1} {1,1,0,0} {0,0,1,1} {1,-1,0,0} 
{1,1,1,-1}{1,1,-1,1}{1,-1,-1,-1}{0,1,0,-1}{1,0,1,0}{1,0,-1,0}{0,1,0,0}{0,0,1,0}{1,0,0,1} 
{1,1,-1,-1}{0,1,0,1}. 

(viii) As shown in [HH|, Peres' 4-dim vectors j22j build a hexagon with 24 vertices 
and 24 edges and not with 22 edges as presented by Tkadlec in [HÜ]- (One can easily 
verify that the edges {1,-1, 1,-1}{1,1,-1,-1}{1,-1,-1,1}{1, 1,1,1} and {1,-1,1,1}{1,-1,-1,-1} 
{1,1,1,-1}{1, 1,-1,1} are missing in the middle of Fig. 1 in |39j-) This Peres' 4-dim 24-24 
KS System contains Systems (a) and (c) from Fig. El and all the Systems from Fig. HJ 

(ix) 4-dim Systems with more than 41 vectors cannot have Solutions from { — 1, 0, 1}, 
and there are no such Solutions to Systems without 0-1 states with minimal loops of 
size 5 up to 41 vectors [there are altogether two such Systems: 38-19(5) given in See. El 
and a 40-20(5) System], what brings the Hasse (Greechie) diagram approach to the KS 
problem [3B1 EU] into question. 

(x) It can easily be shown that a 3-dim system of equations representing diagrams 
containing loops of size 3 and 4 cannot have a real Solution. For loops of size 3, e.g. 
123,345,561 the proof runs as follows. Let us choose 1={1,0,0}, 2={0,1,0}, 3={0,0,1}, 
and i={an, a i2 , a i3 }, i = 4,5,6 and consider block 345. Using 3-5=0 we get a 53 = 0. 
Let us next consider group 561. Using 5-1=0 we get a 5 i = Hence, 5={0,a 52 ,0} and 
is therefore collinear with 2. Thus, the system cannot have a Solution. The proof for 
loops of size 4 is similar, only a little longer. 

(xi) The smallest 3-dim Systems without a 0-1 valuation have a minimal loop of 
size 5, 19 vertices and 13 edges [See. El Fig. El(d)], but they do not have real Solutions. 
We scanned all Systems with up to 30 vectors and 20 orthogonal triads and there are no 
KS vectors among them. This does not mean that Conway-Kochen's system (CK) [26., 
p. 114] is the smallest KS system, though. It turns out that we cannot drop vectors 
that belong to only one edge from orthogonal triads because (a) there are cases where 
a Solution to a füll system allows 0-1 valuation while one to a system with dropped 
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vectors does not and (b) there are cases where the füll System does not allow 0-1 val- 
uation but has no Solution. So, CK is actually not a 31 but a 51 vector system: 123, 
145 ) 267,2AB,3CD,CEF,CGm,DIn J DKL,6EM J 6KN,7I0,7GP,4GQ,4Ko,5Ep,5IS ) ALW,AFX, 
BSY , BQZ , 3cf , 3de , cOh , dMT , cN9 , dP8 , eSl , f Qg , iRl , j kl , iFa , j Lb , kOU , kMV , RPH , RNJ 
with 37 edges. (Tkadlec's claims [HS] that CK can have 55 and 56 vectors and 54 
edges are wrong.) A Solution to CK is 12. . . op={0,0,l}{l,0,0}{0,l,0}{l,-l,0}{l,l,0} 
{0,1,-1} {0,1,1} {2,5,1} {2,5,-1} {0,1,2} {0,2,-1} {1,0,1} {1,0,-1} {1,-1,-1} {1,2,-1} {1,1,-1} 
{2,-l,-5}{l,-l,l}{2,-l,5}{l,l,l}{l,-2,l}{2,l,l}{2,-l,-l}{2,l,-l}{2,-l,l}{l,l,2}{l,2,0} 
{l,-l,-2}{2,-5,l}{2,l,5}{2,l,-5}{5,2,-l}{5,-2,l}{5,l,2}{5,-l,-2}{l,2,5}{l,-2,-5}{l,0,2} 
{l,0,-2}{2,0,l}{2,0,-l}{l,-5,2}{2,-5,-l}{2,-l,0}{2,l,0}{l,-2,0}{l,5,-2}{l,-2,-l}{l,2,l} 
{1,1,-2}{1,-1,2}. Thus, when all the vectors are taken into account, Bub's system [23] 
with 49 vectors and 36 edges: 123 , 345 , 167 , AB6 , AC4 , DEG , DFH , F90 , E8V , 5 JI , 7MN , Gla , 
HNh , 7LT , 5KR , DAe , UTS , PRS , 1 GP , 3HU , 3 V j , Pgh , Uba , 1 i , VZg , OYb , 6Xk , 4Wn , Sde , de i , 
df j , imn , j lk , akQ , hnQ , eQ2 is so far the smallest. 




Figure 5. (a) System with dropped vectors that belong to only one edge (4,6,A,C,D) 
does not have any set of 0-1 states. Taking any of these vectors into account results 
with Systems with at least one set of 0-1 states; (b) Neither the System itself nor the 
Systems obtained by dropping vectors (J or K or both) allow 0-1 states. The latter 
Systems have Solutions while the original System does not have any. 

Let us see why we cannot drop vectors that belong to only one edge in detail. First, 
as mentioned above, if we drop all vectors that belong to only one edge Fig. Ü3(a) we 
get: 123,35, 567, 789, 9B,B1, 28. This system has no 0-1 states. But if we add back 
a Single such vector, say 123,345,567,789,9B,B1 ,28 or 123,35,567,789,9B,B1,2D8, 
the system has at least one set of 0-1 states. All these Systems do have Solutions. The 
opposite Situation is given by Fig. 0(b). The system does not admit 0-1 states but has 
no Solution. If we dropped J or K or both we would have a system with no 0-1 states and 
the Systems would have Solutions. E.g., the system with dropped K has the following 
Solution: 12... IJ = {0,0,0,1}{1,0,0,0}{0,1,1,0}{0,1,-1,0}{1,0,0,-1}{1,-1,-1,1}{1, 1,1,1} 
{1,-1,1,-1}{0,1,0,-1}{1,0,-1,0}{0,1,0,1}{1,-1,1,1}{1,1,1,-1}{1,1,-1,1}{0,0,1,1}{1,-1,0,0} 
{1,1,0,0}{0,0,1,0}{1,-1,-1,0}. Second, in any KS diagram and therefore in Kochen- 
Specker (see Fig. EJ), Bub, Peres, and Conway- Kochen 's ones in particular, only all 
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vectors together make a complete description of their KS sets. Recall that one arrives 
from an n-tuple to an adjoining one by rotation around an (n— 2)-dim subspace. E.g., 
in Fig. |H1 (ii) one starts with 123 and by rotations around 3, 5, and 7 one arrives at 789. 
So, 4 and 6 are indispensable for the construction and cannot be dropped as also shown 
by Larsson [2%] . 




Figure 6. Historical Kochen-Specker 192 (117) graph |23 Lemma 2, p. 68, Fig. p. 69] 
in the MMP diagram notation. Inset (i) shows the hexagon with the adjoining triangle 
from |2U Lemma 1, p. 68]. Inset (ii) shows the same graph in our MMP diagram 
notation: triangles translate as edges and everything eise stays the same except that 
Kochen and Specker drop the vectors that do not share edges, in particular, vectors 4, 
6, A, C and D. In [231 Fig. p. 69] a = p , b = q , and c = r hold. Hcrc we ghie these 
points together graphically. The groups of hexagons p, q, and r here represent the 
hexagons containing pi, and 7"j, i = 0, . . . , 4 in Fig. p. 69]. Inset (iii) represents 
the 27(17)-point graph from |2HI Fig. p. 70] in the MMP diagram notation. 

Special attention is deserved by the first KS graph ever, given by Kochen and 
Specker themselves [2B]- We translated it into the MMP diagram notation in Fig. El 
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where we also explain the correspondence between the two notations. Vectors of the 
KS graph are all contained in the three groups of five hexagons of the type shown in 
the inset (ii) of the figure. Since each such hexagon contains 13 vectors and since two 
hexagons in each group p, q, and r share a vector (vectors a, b, and c, respectively) 
this makes 192 vectors (vertices) and 118 edges. By dropping vectors that do not share 
edges Kochen and Specker obtained 117 vectors. 

Let us just mention here that Kochen and Specker's 27(17)-point graph [Fig-EJ/m/] 
provides a partial Boolean sub-algebra (characterising the Operations of commensurable 
observables in both quantum and classical mechanics) that cannot be embedded into 
a Boolean algebra (i.e., not all classical tautologies from the Boolean algebra corre- 
spond to equalities in the partial Boolean algebra). The graph does allow 0-1 states. It 
has properties similar to those of its hexagons [cf. Fig. EJ-(a)] since it represents a vec- 
tor System 1 23 , 345 , 567 , 789 , 9AB , BC 1 , 4DA , EFG , GHI , I JK , KLM , MNO , OPE , HQN , IRK , 7RE 
with the following components: 12... QR = {{0,1,-2}{5,2,1}{1 -2,-l}{l,0,l}{l,l -1} 
{2 r l,l}{0,l,l}{2,l r l}{l,-l,l}{l,0,-l}{l,2,l}{5,-2,-l}{0,l,0}{0,l r l}{2,l,l}{l r l r l} 
{1,1,0} {1,-1,2} {5,1,-2} {0,2,1} {5,-1,2} {1,1,-2} {1,-1,0} {1,1,1} {2,-1,-1} {0,0,1} {1,0,0}} 
which does have a set of 0-1 states when all vectors are taken into account and does 
not have it when the vectors that do not share edges are dropped. However, we obvi- 
ously cannot dispense with vectors that build the System and therefore we cannot use 
this System for proving the KS theorem. Of course, Kochen and Specker were aware of 
this fact too and this is why they designed the aforementioned 15-hexagon 192-vector 
System. 

(xii) The concept of KS dual diagrams plül l4"ül l4"Tj is apparently either a misnomer 
or insufficiently defined. Tkadlec claims that one arrives from a Standard 3-dim 
KS diagram to its dual so as to "replace the role of points [vertices] and smooth 
curves [edges]: points [vertices of the dual diagram] represent blocks [edges of the 
Standard diagram] and maximal smooth blocks [maximal (?) edges of the dual diagram] 
represent atoms [vertices of the Standard diagram]." [3Ü] The Instructions for such a 
construction are ambiguous, but two figures are given in IHH] and and we have 
tested them. One is a dual diagram to Peres' 57 (33) diagram % and it reads: 
123 , 345 , 567 , 869 , 9 AH , 8C2 , 7DG , HG1 , 4BA , CBD , 6gE , BhE , 31 J , 2R0 , 1VU , VPN , UML , JKN , 
0KL , IQM , RQP , j SK , j iQ , UWX , Veb , Gf a , HCZ , ZYb , XYa , WdC , edf , TFd , TcY , IkE , 11 j , lmT 
and the other: 123 , 145 , 16C , 768 , 7HK , 4FB , GEC , 89A , 51 J , HGI , EF9 , KBD , JAD , CDV , KLM , 
BON , DgS , VUT , SP2 , QRS , MQU , NPT , c6R , GPd , VWX , X3Y , 3Ze , EZQ , ab J , Yf A , cde , ehD , acW 
is dual to Conway-Kochen's diagram. Of these two diagrams only the latter does not 
allow 0-1 state. The former has at least one set of 0-1 states. Then our solvers prove 
that the equation system corresponding to the latter diagram does not have a Solution. 
Hence, neither of the two diagrams is a KS set, and we would expect a KS dual to be a 
KS set. 

(xiii) We obtain the class of all remaining (non-KS) vectors from 7i n by first filtering 

f 123,39R,89A,47D,56E,DRE,EFG,CBD,NML,LKE,Djq,QST,PJI ) HK0,RVX,RUW,14Y,lZ5,4aA,5b8,8gB J 
AhF.TcH.edl.CiD.GjP.reM.efS.ClN.GkT.NqX.PsV.DrU.MmU.SnV.HoX.IpW.TtW^uB^vF 
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the MMP diagrams so as to keep only those that allow 0-1 states. Out of these, a second 
filter then keeps only those diagrams whose corresponding equations have Solutions. 
Vectors corresponding to these Solutions are the wanted non-KS vectors. 

(xiv) The presented algorithms can easily be generalised beyond the KS theorem. 
One can use MMP diagrams to generate Hilbert lattice counterexamples, partial Boolean 
algebras, and general quantum algebras which could eventually serve as an algebra for 
quantum Computers. jlHj One can also treat any condition imposed upon inner products 
in W 1 to find Solutions not by directly solving all nonlinear equations but also by first 
filtering the corresponding diagrams and solving only those equations that pass the 
Alters. 
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